Triplet-singlet relaxation in semiconductor single and double quantum dots 
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We study the triplet-singlet relaxation in two-electron semiconductor quantum dots. Both single 
dots and vertically coupled double dots are discussed. In our work, the electron-electron Coulomb 
interaction, which plays an important role in the electronic structure, is included. The spin mixing 
is caused by spin-orbit coupling which is the key to the triplet-singlet relaxation. We show that the 
selection rule widely used in the literature is incorrect unless near the crossing/anticrossing point 
in single quantum dots. The triplet/singlet relaxation in double quantum dots can be markedly 
changed by varying barrier height, inter-dot distance, external magnetic field and dot size. 
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I. INTRODUCTION 

The application of semiconductor quantum dots (QDs) 
in generating spin-based qubits^ is one of the focuses in 
the field of spintronics.— There are two types of qubits 
investigated extensively recently— One is based on the 
transition between single-electron Zeeman sublevels 4 ^ 5 - 
and the other is based on two-electron triplet-singlet 
(TS) st at es 1 1 1 1 Among these works, the decoher- 
ence time of the spin states, including both the spin de- 
phasing time 7 -^ and spin relaxation time; 5 -' 9 - ' 10 ' 11 has at- 
tracted much attention as a thorough understanding of 
it is one of the prerequisites of the application. There 
are many works on spin relaxation reported, especially 
in single-electron QD s 3>M^M^ffi^ Re- 
cently, the TS relaxation time of two-electron system 
has also been investigated ! 24 ' 25 ' 26 ' 27 It was proposed that 
various mechanisms, such as the electron-phonon scat- 
tering together with the spin-orbit coupling ) 28 ' 29 the hy- 
perfine interaction , 30 ' 31 and the cotunneling effect, could 
induce TS relaxation^ However, the mechanism involv- 
ing electron-phonon scattering is usually treated as the 
key one because the nuclei-mediated relaxation^ 2 , and the 
cotunneling can be weakened via tuning external mag- 
netic field and tunneling rates?^ respectively^ 4 - Specif- 
ically, Climente et al. used exact diagonalization tech- 
nique to calculate the two-electron spectrum structure 
and the phonon induced TS relaxation in parabolic sin- 
gle QDs^ 4 - They demonstrated the crucial role of the ex- 
cited states on spectrum structure and showed a slow 
decrease of the relaxation time away from the TS cross- 
ing in contrast to a sharp increase in the vicinity of the 
crossing point, when the magnetic field is increased from 
zero Tesla. This feature agrees qualitatively with the re- 
cent measurement i 1 - 1 - Furthermore, their results indicated 
that the spin-down triplet state coupled with the singlet 
ground state through the spin-orbit coupling has a much 
shorter lifetime compared to the other two triplet states. 
This was understood by the so called "selection rule" 
based on the perturbation using the lowest two single 
electron levels. Similar perturbative discussion was also 



given in Ref. [25|. Meunier et al. obtained perturbative 
wave functions from the selection rule and treated the 
spin-orbit coupling coefficient as a fitting parameter iii 
Using these functions, they fit their experiment data 
with clcctron-phonon-scattering-induced TS relaxation 
and obtained a particularly small spin-orbit coupling co- 
efficient. They attributed the reduction of the coupling 
coefficient to the neglect of high excited states. Sasaki et 
al. pointed out that the selection rule was correct only in 
the vicinity of the TS crossing point ^ which seems to be 
more correct intuitively. According to the previous work 
by one of the authors 1 ^ and confirmed by Destefani and 
Ulloaji 7 - the spin-orbit coupling in quantum dots is very 
strong and a large number of basis functions are needed 
in order to achieve convergence even for the lowest few 
states. Therefore, whether the selection rule based on the 
lowest few levels remains unchanged when many upper 
levels are involved remains questionable to us. There- 
fore, in this work we will first reinvestigate the selection 
rule based on exact diagonalization method, jointly with 
perturbation method with many basis functions. 

The investigation on TS relaxation in double QD ar- 
chitectures is very limited. Recently, Wang and Wu stud- 
ied the single-electron spin relaxation in vertically cou- 
pled double QDs and showed that the spin relaxation can 
be efficiently manipulated electronically by the inter-dot 
barrier.— This suggests that the two-electron TS relax- 
ation should also be manipulated by tuning inter-dot bar- 
rier height. This is another issue we are going to explore 
in this work. 

We organize the paper as following. In Sec. [IT] we set 
up the model and lay out the formalism. Then in Sec. lIIII 
we show our numerical results. We discuss the single dot 
case in Sec. MI Al We first show the exact diagonalization 
results with sufficient basis functions. We then reexamine 
the selection rule by using more basis functions instead 
of the lowest two, both perturbatively and exactly. We 
show the selection rule widely used in the literature is not 
correct except near the TS crossing/anticrossing points. 
In Sec. IIII B[ we show the results of double QDs. We 
summarize in Sec. IIVI 
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II. MODEL AND FORMALISM 



and 



We start our investigation from a vertically-coupled 
double QD: Electrons are confined by a parabolic po- 
tential V c (x, y) = \m*uj 2 {x 2 + y 2 ) (corresponding to the 
effective dot diameter rf = \/hir/m*LOo) along the x- 
y plane ; 33 ' 34 with m* representing the effective mass. 
Along the z-axis, a strong confinement is given by 




\z\ < \a , 

\a < \z\ < |a + d , 
otherwise , 



(1) 



with Vo, the inter-dot barrier ^ By taking a = 0, one 
comes to the single dot configuration. The single-electron 
Hamiltonian with magnetic field along the growth direc- 
tion (z) is given by 



2m* 



V(r)+H so (P)+H 2 



(2) 



in which F(r) = V z (z) + V c (x, y) and P = -iHV + e/cA 
with A = (B/2)(— y, x, 0). H so represents the spin-orbit 
coupling which is the key to the spin flip. In this work, we 
only consider the Dresselhaus spin-orbit coupling^ as the 
Rashba coupling^ is comparably small in GaAs QDs. 36 
Hence H so = 7 h a, with h = [P X {P 2 - P 2 ),P y {P 2 - 
P 2 ), P Z (P 2 - P 2 )}£- For small well width, it reduces to 

H so = ^(P 2 )(-P x a x +P y a y ) , (3) 

with (P 2 ) the average of P 2 over the electronic states 
defined by V z (z). Hz — \g^BBcr z is the Zeeman split- 
ting with g being the Lande factor. We define H = 
+ V(r), whose eigenvalues and eigenfunctions can be 
obtained from the Schrodinger equation 

Ho\<t>) =E\<t>) . (4) 

Previous work on single-electron QDs gives the solution 
of the lateral part of this equation ) 14 ' 33 ' 34 where the exact 
energy levels are given by 

(5) 

The wave 



E nl = m(2n + \l\ + 1 ) + hluj B 



with = y/uj 2 + ujg and ojb = eB/(2m* 
functions read 

(r\nl) = N n ^(ar)^e~ {ar)2/2 L^((ar) 2 )e ue , (6) 

with N n ,i = (a 2 n\/ir(n+ l^l)!) 1 / 2 and a = ^/m*n/ti. 
is the generalized Laguerre polynomial. In these equa- 
tions, n = 0,1,2,... is the radial quantum number and 
/ = 0, ±1,±2, ... is the azimuthal angular momentum 
quantum number. By solving the z-component of Eq. 
Q, we obtain the lowest two electronic states along the 
z-axis as following: 



C$sm[k(z- f -d)\ , 
C 2 ° cosh(/3z) , 
C%sm[k{-z- | -d)] , 



\z\ < ~a , 



d < z < 



C\ sin[fc(z - 
G\ sinh(/3z) 
C\ sin[fc(z -\ 



f-d)] 



< z < 
I £ ha 



(8) 



§ + d)], -§-d<z<-§ 



in which k 2 = 2m*E z /h 2 and (3 2 = 2m*(V - E z )/h 2 
with E z denoting the energy along this direction. We 
use the superscripts "0" and "I" to denote the even and 
odd parity respectively. The total spatial wave function 
is then denoted by \nln z ), with n z = and I in this work 
to distinguish the above even and odd states along the 
z-axis. Due to the strong confinement along the z-axis, 
levels higher than n z = 1 are neglected. It is noted that 
when we refer to the single QDs, we only keep the lowest 
state (the even one) due to the small well width. 

For two-electron system, the total Hamiltonian is writ- 
ten as 



Htot — {H e + H 2 + He) + H\ v + H ep + H p 



(9) 



de- 



In this equation, the third term Hn = i f r 

^ ' ° 47re K;|ri— r 2 | 

scribes the Coulomb interaction between the two elec- 
trons with k representing the static dielectric constant. 
H p = Tiqxhuiqxa^aqx represents the phonon Hamilto- 
nian, and H ep = Sq^Mq^a^ + a q A)exp(iq ■ r) is the 
Hamiltonian of the electron-phonon interaction. The su- 
perscripts "1" and "2" label the two electrons. 

We construct two-electron basis functions from the sin- 
gle electron wave functions. To see the physics clearly, we 
construct our two-electron basis functions in either sin- 
glet or triplet forms. Taking two single-electron spatial 
wave functions \nilin z i) and \n 2 l 2 n z2 ) (denoted as \Ni) 
and \N 2 ) for short) as an example, the singlet functions 
can be constructed by 



\s) = (I n> - 1 it)) 



Ni = N 2 , 



(10) 



(7) 



\\{\NxN 2 ) + \N 2 N X )) , N X ^N 2 , 
and the triplet functions for N x ^ N 2 by 

\T+) = -^m^) - \N 3 Ni)) ® | TT> , (11) 
|T ) = ~(|tfiJVa) - \N 2 N x )) ® (| U) + I IT)) ,(12) 
\T-) = -^mNz) - \N 2 N X )) ® | U) . (13) 

Here, N and N', in the ket \NN'), represent the spa- 
tial quantum numbers of the first and the second elec- 
trons respectively. We define the total angular momen- 
tum L = li + l 2 and denote the total spin (S, S z ) with 5*2 
representing the z-component of the total spin S. 

Then, we calculate the matrix elements of the Coulomb 
interaction and the spin-orbit coupling^ under these ba- 
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sis functions. The Coulomb matrix elements can be ex- 
pressed in the form 



°In 1 +In 2 ,I n i +l N * 



X Q(N-l,N 2 ,N' v N2) 



(14) 



in which Q is given in detail in Appendix [XJ Thus we 
obtain the two-electron Hamiltonian. By diagonalizing 
the two-electron Hamiltonian, one obtains all the en- 
ergy levels and eigenfuctions. We identify a state as sin- 
glet/triplet if its amplitude of singlet/triplet components 
is larger than 50 %. We rewrite the spin-orbit coupling 
Hamiltonian [Eq. ([3])] using the ladder operators asS^ 



H so = lc (P+S+ + P~S- 



(15) 



with the coupling coefficient 7 C = Jj(P^). Then it is 
noted that P and 5 change L and S z by one unit, 
respectively. It suggests that a state with (L,S Z ) can 
only be coupled with the states with (L + 1, S z + 1) and 
(L- 1,5,-1). 

Treating \i) and |/) as the initial and final states, we 
can calculate the phonon-induced relaxation rate from 
the Fermi Golden Rule 



r W = T El A Vl 2 |</lxN)| 2 [^A<5( e/ - ei -^ qA 

qA 
+ (™qA 



l)S(e 



f 



qAj 



(16) 



in which x(q,ri,r2) = e' qri + e lqr2 comes from the 
total electron-phonon interaction Hamiltonian H ep — 
H\ p + H 2 . Here, n q A represents the Bose distribution of 
phonon with mode A and momentum q. In our calcula- 
tion, the temperature is fixed at K. Therefore only the 
phonon emission process occurs. 



III. NUMERICAL RESULTS 



calculation, we include 
phonon scattering due 



In the numerical 
the electron-acoustic 
to the deformation potential with |Af qs ;| 2 = 
HE 2 q/2Dv si r^and due to the piezoelectric field 



"si! 

with |M qp/ | 2 = (32hTT 2 e 2 el A /K 2 Dv s i)[(3q x q y q z 
for the longitudinal mode^ 9 . and £j=i,2|Afqpt- 1 2 = 

{i2hn 2 e 2 e 2 jK 2 Dv st q 5 )[q 2 x q 2 v W y qlW z qlAM x qyqz)VQ 1 \ 
for the two transverse modest Here, S = 7 eV stands 
for the acoustic deformation potential; D = 5.3 x 10 3 
kg/m 3 is the GaAs volume density; ei4 = 1.41 x 10 9 
V/m denotes the piezoelectric constant and the static 
dielectric constant k is 12.9; v s i = 5.29 x 10 3 m/s 
corresponds to the longitudinal sound velocity and 
v s t = 2.48 x 10 3 m/s corresponds to the transverse 

41 

onei^i 

In our calculation, g factor is — 0. 11 i and the Dres- 
selhaus coefficient 7 is 21.5 A 3 -eV42 The typical electron 
effective mass m* in GaAs is 0.067mo^ with mo being 
the free electron mass. 
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FIG. 1: (Color online) (a) The lowest four energy levels vs. 
magnetic field B in single QD. The TS anticrossing point 
between T_ and S is shown and the range near this point is 
enlarged in the inset, (b) r -1 of the three transition channels 
vs. the magnetic field. In the calculation, d =5 nm and 
do = 30 nm. 



A. Single dot 

We first set a = to investigate the single dot case 
by exact diagonalization method with the lowest 800 sin- 
glet and 2220 triplet basis functions. Under the basis, 
the energy levels and the TS relaxation rates are well 
converged. The magnetic field dependence of the first 
four levels and that of the TS relaxation rates are plot- 
ted in Fig. [T] In the calculation, we take the well width 
2c? = 10 nm and the effective diameter do — 30 nm. 
From Fig. [T{a), one notices that the ground state is a 
singlet denoted as "S" , in a wide range of the magnetic 
field (from T to 2.6 T approximately). In this region, 
the first three excited states are triplet states, labeled 
as |T + ) (spin-up), |T ) (spin-zero) and |T_) (spin-down), 
and the energy of |T_) is the highest one among the three 
because of the Zeeman effect. When the magnetic field 
increases from 2.6 T, one further observes a TS crossing 
between the singlet and the two triplets (|T+) and \Tq)). 
Moreover, a TS anticrossing point (with a small energy 
gap shown in Fig. [Ha)) a l so exists between the singlet 
and |T_) triplet state due to the Dresselhaus spin-orbit 
coupling. From the calculation, we notice that the ma- 
jor components of \S), |T±) and |T ) are |5 X ), \T£) and 
ITq 1 ), which are the lowest singlet and triplet basis func- 
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tions. Specifically, using the lowest two single-particle 
wave functions, \nln z ) with n = n z = 0, / = and — 1, 
one can construct IS 1 ) with |000) and |000), and |T^) and 
IT 1 ) with 1 000) and |0 - 10) according to Eqs. (flt ty -ip ]) . 
Therefore, the quantum numbers (L,S Z ) of IS 1 ), |Tf), 
ITq 1 ) and IT 1 ) are different, i.e., (0,0), (-1,1), (-1,0) 
and (—1,-1) respectively. 

From Fig. [2(b) one observes that the TS relaxation 
rates increase slowly with the magnetic field away from 
the crossing/anticrossing points, but decrease dramati- 
cally in the vicinity of the crossing/anticrossing points, 
in agreement with the measurement qualitatively* 1 ^ The 
relaxation rate reaches maximum where the wave length 
of the emissive phonon is comparable with the dot size^ 3 - 
In our calculation, the TS splitting, i.e., the energy be- 
tween the triplet and the singlet, Ats ~ 0.2 meV. The 
corresponding half-wavelength of the transverse phonon 
is therefore about 30 nm as the dot diameter do. This 
feature was interpreted as the competing effects of the 
magnetic field on the elcctron-phonon coupling and the 
spin-orbit coupling, 2 = Actually, the strength of the spin- 
orbit coupling is proportional to a [see Eq. (j6))] which 
increases with the magnetic field ^ whereas the electron- 
phonon scattering becomes rather weak when the emis- 
sive phonon momentum decreases^ 

Surprisingly, our results are very different from those 
shown in the previous work, where the transition rate of 
|T_) is much larger than those of the other two triplet 
states \T + ) and \Tg)^ In that work, the authors in- 
terpreted their results by the selection rule based on 
the perturbation method including the lowest four ba- 
sis functions, i.e., IS 1 ), |T|), IT,, 1 ) and IT 1 ). Under that 
basis, only |T_) is coupled with \S) through the Dres- 
selhaus spin-orbit coupling according to Eq. (fTSJ). So 
only the transition from |T_) to \S) can occur. Thus 
they concluded that the transition rate from |T_) to \S) 
is much larger than those of the other channels even 
though much more (instead of four) basis functions are 
included. In fact, this selection rule is widely used in the 
literaturei 10 ' 11 However, as one needs many basis func- 
tions to achieve convergence even in the single-electron 
QD system^ whether the selection rule from the low- 
est four basis functions is robust against the inclusion of 
higher basis functions remains an open question. Here 
we reexamine the selection rule with more basis func- 
tions. Assuming the perturbation based on the lowest 
four states (S 1 ), \T±) and \Tq) is adequate to describe the 
real physics, we expect the selection rule should always 
be valid when more basis functions are included. Specif- 
ically, we now use four single-electron functions |000), 
|0 - 10), |010) and |0 - 20) to construct the two-electron 
basis functions. Keeping only the index of I from \nln z ) 
since the other two are fixed, the six lowest singlet states 
are constructed by |0)|0), |0)| - 1), |0)|1), | - 1)| - 1) , 
1 0)| — 2), and | — 1)|1) separately and the three lowest 
triplet states are constructed by |0)| — 1) in the way of 
Eqs. (fTU)) - (fTB|) . We denote these nine basis functions as 
IS 1 ), |S 2 ), IS 3 ), |S 4 ), |S 5 ), |S 6 ), {T 1 }, IT 1 ) and |T±) in 



sequence, and the quantum numbers (L,S Z ) are (0,0), 
(-1,0), (1,0), (-2,0), (-2,0), (0,0), (-1,1), (-1,0), 
and (—1,-1) respectively. Therefore, only the singlet 
states IS 1 ) and |S 6 ) can mix with |Tl); |S 4 ) and |S 5 ) 
can mix with \T+), according to Eq. <|T~5|> under these 
basis functions. No mixing occurs to the state ITq 1 ). 
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FIG. 2: (Color online) Comparison of the TS relaxation rates 
of T- and T+ in single QD calculated from the perturba- 
tion method and the exact diagonalization method limited 
within nine basis functions (IS" 1 ) to |S 6 ), |7±) and \T^)). 
Dashed curve and A: 1/tt + ->s from the exact diagonaliza- 
tion and perturbation methods respectively; Solid curve and 
• : 1/tt_^s from the exact diagonalization and perturbation 
methods. 



As the Coulomb interaction is too strong to treat per- 
turbatively, we first diagonalize the Hamiltonian with the 
Coulomb interaction included to obtain a new set of basis 
functions, i.e., IS 1 ) = a^S 1 ) + 61JS 6 ), |S 2 ) = ^S 1 ) + 
6 2 |S 6 ), |S 3 ) = a 3 |S 4 ) + 6 3 |S 5 ), |S 4 ) = a 4 |S 4 ) + 6 4 |S 5 ), 
|S 5 ) = |S 2 ), |S 6 ) = |S 3 ), |?|) = |T|), IT 1 ) = IT 1 ), and 
|Ti) = {T 1 ). Here a* and bi are obtained from the nu- 
merical diagonalization. The corresponding eigenvalues 
are E\ to Eq, E + , Eq and E-, respectively. Then we 
treat the spin-orbit coupling as perturbation under the 
new basis functions. The lowest four states then read 



\To) 
IT- 
IS) 



> = in) 

= ITq 1 ) 

> = |t±) 

= 0l\S 1 ) 



d\\S 4 ) 

elis 1 ) 

-0 2 |S 6 ) 



? 2 _|S 5 ) 



■0 2 |S 6 ) , 
■&.\Tl) , 



(17) 
(18) 
(19) 
(20) 



with 



9f 

a'A+b'B 



_^ (b* -<) A p,l(2) 



e 



1(2) 



'<+- 

and (9 3 



jl(2) _ @ 1(2) q1(2) _ 



aiA+biB 
E--Ei ' 



Here 9^ 
-^ 7 *a(l + 



ai(bi), A — — ij*a(l ~ eB/2ha 2 ) and B ■ 

eB/2ha 2 ) with 7* being j(Tr/2d) 2 . 

Obviously, the transitions from both \T + ) and \T_) to 
|S) can occur according to Eqs. (|17 p - (|2"0)) . The matrix 
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elements |(/|xN)| 2 m Eq. (| 16[) now read 

|(5|x|T+)| 2 = |flJX<5 1 |x|S 6 ) +ef^<5 6 |x^ 4 )| 2 

= (zi) 2 |£| 2 / 2 fe), (21) 
|(S|x|T_)| 2 = lO-^lxIS 1 ) + *f£ <S 6 |x|S 6 ) 

+0f (rl| x |rl)| 2 / 2 ( gz ) 

= |2tCi-teC 2 | 2 / 2 fe) , (22) 
with ;r = fc 2 /4a 2 , i = e _a: , I(q z ) — tt 2 sm(dq z )/{dq z [n 2 — 

(d qz ) 2 }}, a = 0f + C^i + efel , c 2 - of + 2e 2 *e 2 _ 

and £ = ^s*^ 2 + \[2Q 2 * d\. We calculate the relaxation 
rates of these two channels and plot the results in Fig. 
[21 One notices that the two sets of dots (blue ▲ for \T + ) 
and red • for |T_)) are quite close to each other and even 
show a crossing. In other words, the selection rule is vi- 
olated. We also present the exact diagonalization results 
under the same basis functions IS 11 )-^ 6 ), \Tq) and |T±) 
in Fig. [5] (blue dashed curve for \T + ) and red solid curve 
for |T_)). It is seen that the diagonalization results al- 
most exactly match the perturbation results. This match 
further confirms that both our exact diagonalization and 
the perturbation calculations are correct. Compare Fig. 
[2] with Fig. |Tfb) , it is obvious that the high excited lev- 
els manifest themselves markedly in the relaxation rates. 
From our calculation, we notice that the coefficients in 
Eqs. (fTT)) and ([TO)) are comparable. This is because that 

1(2) 

the denominators E- — Ei in 9_ are close to E+ — E{ 

in 0+ ' . This explains the reason why the curve of |T_) 
is close to that of \T + ) in Fig. 

However, it is noted that the selection rule works well 
in the vicinity of the crossing/anticrossing points both 
in Figs.[]Jb) and This can be understood from Eqs. 
([T7|) and (fT9|) . Near the TS crossing point where E + ~ 
E\, the energy splitting E + — E^a is finite. Therefore 

1(2) 

# + only changes slightly compared with the region away 
from the TS crossing. Similar is true for the coefficients 

1 ( 2) 

of \Tq). In contrast, 8_ ' is very large when E- ~ E\. 
Therefore the transition rate from |T_) would be much 
larger than those from \T+) and \Tq), i.e., the selection 
rule is valid in the vicinity of the TS crossing/anticrossing 
point. Moreover, the effect of the Zeeman splitting also 
makes the transition rate of |T_) larger than those of 
|T_i_) and \Tq) because of the larger phonon momentum q. 
Specifically, the energy splitting between |T_) and \S) is 
about 0.18 meV at B — 2.5 T in Fig.rjja), which is much 
larger than that between \T+) (|T )) and \S), i.e., ~ 0.06 
meV (0.12 meV). As the transition rates are proportional 
to q m with to > varying for different mechanisms, the 
rate of |T_) is much larger than those of \T + ) and \Tq). 

B. Double dot 

Now we turn to study the TS relaxation rate in weakly 
coupled double QDs using a basis functions including 400 
singlet and 1080 triplet states. In the calculation, a = 8 



nm and d = 7 nm. In this part we still use \T±) and \Tq) 
(| S)) to denote eigenfunctions of the lowest three triplet 
states (lowest singlet state). To determine the contri- 
bution of the energy levels along the z-axis, we take the 
barrier height Vq = 0.25 meV, the lowest one in our calcu- 
lation, as an example. In this configuration the splitting 
between the first and the second levels along the z-axis is 
about 1 meV and that between the second and the third 
levels is much larger, about 0.2 eV. Compared with the 
lateral confinement (~ 4 meV for do = 30 nm), we only 
need to include the lowest two in our calculation. 
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FIG. 3: (Color online) r 1 vs. barrier height in double QDs. 
In the calculation, a — 8 nm, d — 7 nm, do = 30 nm and 
B = 0.5 T. 

We first investigate the TS relaxation rate as a func- 
tion of the barrier height. In the calculation, d = 30 
nm and B = 0.5 T. As shown in Fig. [31 each transition 
rate first increases slowly until it reaches the maximum 
around V ~ 0.28 eV where the TS splitting A T s ~ 0.4 
meV corresponding to the wavelength of the emissive 
phonon being comparable with the dot size do^ After 
that, the TS relaxation rate decreases rapidly with the 
barrier height. This would offer us a scheme to manip- 
ulate the TS relaxation in double QDs. Similar features 
(not shown here) are obtained when we increase the inter- 
dot distance. The dramatic decrease of the relaxation 
rate can be understood as following. When the barrier 
height becomes higher or the inter-dot distance becomes 
larger, the inter-dot coupling is weakened and the energy 
splitting between the lowest two levels along the z-axis 
becomes smaller. As a result, the splitting between \T±) 
(| To)) and \S) decreases too. This causes the decrease 
of the TS relaxation rate as discussed in the previous 
subsection . 

To have a look at the role of the magnetic field, we 
calculate the average relaxation rate 1/f = (1/t<t + ^s + 
1/tt_->s + 1/ t t ^s)/3 as function of the barrier height at 
different magnetic field in Fig.Q](a), but with the dot size 
do = 30 nm fixed. It is seen from the figure that higher 
magnetic field leads to relatively larger transition rate. It 
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V (V) d (nm) 



TS relaxation rates based on nine basis functions. Com- 
parable transition rates of |T+) and |T_) are immediately 
obtained away from the TS crossing point. The pertur- 
bation results are in good agreement with the exact di- 
agonalization results under the same basis functions. We 
also show that the selection rule works well in the vicinity 
of the TS crossing/anticrossing point due to the effects 
from the Zeeman splitting and the anticrossing. For the 
double QD case, we demonstrate that the TS relaxation 
rates vary more than two orders of magnitude by tun- 
ing the inter-dot barrier. This offers a feasible scheme to 
manipulate the TS relaxation in double QDs. The relax- 
ation rates also sensitively depend on the dot size and 
magnetic field. 



FIG. 4: (Color online) (a) Average TS relaxation rate f _1 vs. 
barrier height at different magnetic fields in double QDs. In 
the calculation, a = 8 nm, d = 7 nm and do = 30 nm. (b) 
r _1 vs. effective dot diameter. In the calculation, a — 8 nm, 
d = 7 nm, V = 0.35 V and B = 1 T. 



is due to the enhanced spin-orbit coupling in strong mag- 
netic field. The influence of the effective diameter of QDs 
with V = 0.35 V and B = 1 T is also shown in Fig. GOT)). 
One finds the transition rates increase with the effective 
diameter do- The reason lies on the different symmetry 
properties of the singlet and triplet states. For the sin- 
glet state, the inter-electron distance decreases with the 
decrease of the dot size. The coulomb repulsion there- 
fore lifts the corresponding energy levels. However, the 
energy lifts of the triplet states are smaller due to the 
antisymmetry property of the triplet states which pre- 
vents the electrons to be close to each other. Therefore, 
the TS splitting becomes smaller with the decrease of 
the dot size. This leads to the rapid decrease of the TS 
relaxation rates. 



IV. SUMMARY 

In summary, we have investigated the TS relaxation in 
single and double QDs. For the single dot case, we find 
that the average relaxation rate first slowly increases with 
magnetic field until it reaches the maximum where the 
wavelength of emissive phonon is comparable with the 
dot size. Then it drops sharply. This result qualitatively 
agrees with the recent measurement ^ Furthermore, our 
result shows the transition rates of the triplet \T + ) and 
\Tq) can be comparable with that of |T_), which violates 
the selection rule in the literature^ We show that the 
selection rule obtained from the lowest four basis func- 
tions does not hold in general cases where much more 
basis functions are needed to converge the triplet /singlet 
states. This is shown perturbatively by calculating the 
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APPENDIX A: Q IN COULOMB INTERACTION 



Following Ref. |1J, we obtain Q in Eq. (|T4|) as 



Q(N 1 ,N 2 ,N[,N^) = j dkfyP NuN > (fey )P K , N2 (fey) 

, A1 , 



where Pn.n' and Wn,n' come from the lat- 
eral and vertical parts of the matrix element 
(n, /, n z \ exp(ik • r)|n', I', n' z ), respectively. P is given 



PN,N>(h) = 



i\n'\ 



y (n+|Z|)!(n' + |r|)! eXP V 4a 2 
" l l' ( ^ 2 \ 



i=0 j=0 



Aa 1 



x(sgn(Z'-Z)^ 



(A2) 



with (7^ = andn = i+j+ (\l\ + \l'\-\l'-l\)/2. 

sgn(a;) represents the sign function. W reads 



Wn,n> = (n z \exp(ik z z)\n' z ) 



(A3) 
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